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Abstract 

, We consider the offset-deconstruction problem: Given a polygonal shape Q with n vertices, can it be 

^SJ ' expressed, up to a tolerance e in Hausdorff distance, as the Minkowski sum of another polygonal shape 

P with a disk of fixed radius? If it does, we also seek a preferably simple-looking solution P; then, 
^ ', P's offset constitutes an accurate, vertex-reduced, and smoothened approximation of Q. We give an 

' 0{n log n)-time exact decision algorithm that handles any polygonal shape, assuming the real-RAM model 

of computation. A variant of the algorithm, which we have implemented using the CGAL library, is based 
on rational arithmetic and answers the same deconstruction problem up to an uncertainty parameter S; 
its running time additionally depends on S. If the input shape is found to be approximable, this algorithm 
, also computes an approximate solution for the problem. It also allows us to solve parameter-optimization 

problems induced by the offset-deconstruction problem. For convex shapes, the complexity of the exact 
decision algorithm drops to 0(n), which is also the time required to compute a solution P with at most 
C/2 , one more vertex than a vertex-minimal one. 
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1 Introduction 

>. 

00 , The r-ojfset of a polygon, for a real parameter r > 0, is the set of points at distance at most r away from 

■ the polygon. Computing the offset of a polygon is a fundamental operation. The offset operation is, for 

instance, used to define a tolerance zone around the given polygon [11] or to dilute details for clarity of 
graphic exposition [161 1201 14] . Technically, it is usually computed as the Minkowski sum of the polygon and 
a disk of radius r. The resulting shape is bounded by straight-line segments and circular arcs. However, 
a customary practice is to model the disk in the Minkowski sum with a (tight) polygon, which yields a 
piecewise-linear approximation of the offset. Our study is motivated by two applications, where such an 
approximation forms the legacy data which a program has to deal with — the original shape before offsetting 
is unknown. This leads us to the question what is the original polygon whose approximate offset we have at 
' hand. Of course, finding the exact original polygon, or even its topology, is impossible in general, because 

^ I the offset might have blurred small features like holes or dents. However, a reasonable choice can lead to a 

more compact and smooth representation of the approximate offset. 

The first relevant problem concerns cutting polygonal parts out of wood. A wood-cutting machine, which 
can smoothly cut along straight line segments and circular arcs, is given a plan to cut out a certain shape. 
This shape was designed as a polygon expanded by a small offset, but with circular arcs approximated by 
polygonal lines comprising many tiny line segments. Thus instead of moving smoothly along circular arcs, 
the cutting tool has to move along a sequence of very short segments, and make a small turn between every 
pair of segments. The process becomes very slow, the tool heats up, and occasionally it causes the wood to 
burn. Moving the cutting tool smoothly and fast enough is the way to keep it cool. If this were the only 
issue, other smoothing techniques like arc-spline approximation [Sj [12] may have been applicable. However, 
we may also wish to reduce the offset radius if a more accurate cutting machine is available — in this case, it 
seems desirable to find the original shape first and then to re-offset with a smaller radius. 



*A preliminary and partial version appeared in the Proceedings of the 27th annual ACM Symposium on Computational 
Geometry (SoCG 2011) 
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A motivation to study this question from a different domain is to recover shapes sketched by a user of a 
digital pen and tablet. The pen has a relatively wide tip, and the input obtained is in fact an approximate 
offset (with the radius of the pen tip) of the intended shape. The goal is to give a good polygonal approxima- 
tion of the intended shape. More broadly, as the offset operation is so commonplace, it seems natural to ask, 
given only an (approximated) offset shape, what could be the original shape before the offsetting. Therefore, 
we pose the (offset- )deconstruction problem which comes in two variants: 
Problem 1: the decision problem 

Given a polygonal shape Q, and two real parameters r,e > 0, decide if there exists a polygonal shape 
P such that Q is within (symmetric) Hausdorff-distance e to the r-offset (i.e., offset with radius r) of 
P 

Problem 2: finding a solution 

If the answer to Problem 1 is YES, compute a polygonal shape P with the desired property. We refer 
to P as a solution of the deconstruction problem. Note that P might be disconnected, even if Q is 
connected (Figure [LT|) . 

Problem 1 can be seen as a special case of the Minkowski decomposition problem which asks whether 
a set can be composed in a non-trivial way as the Minkowski sum of two sets — disallowing a summand to 
be a homothetic copy of the input set. A general criterion for decomposability of convex sets in arbitrary 
dimension has been presented in |19] . A particularly well-studied case are planar lattice polygons, because 
of their close relation to problems in algebra, for instance, polynomial factorization [17]. It has been shown 
that deciding decomposability is NP-complete for lattice polygons . In [B] , decomposability is investigated 
under the constraint that one of the summands is a line segment, a triangle, or a quadrangle. However, 
all these approaches discuss the exact decomposition problem; our scenario of being Hausdorff-close to a 
particular decomposition seems to not have been addressed in the literature. Allowing tolerance raises 
interesting algorithmic questions and at the same time makes the tools that we develop more readily suitable 
for applications, which typically have to deal with inaccuracies in measuring and modeling. 



Figure 1.1: For a given Q, the red P is a candidate summand whose exact r-offset is shaded. Left: For a 
given £, deconstruction is ensured iff 0i < e and 02 < £• Note that, when r decreases, (f)i decreases, but 02 
increases. Middle: Example where Q can be approximated by an r-offset of a P that has much fewer vertices 
than Q. Right: Example where Q can be approximated by the r-offset of a disconnected shape P. 

Contributions. We first present an efficient algorithm to decide Problem 1: For a shape Q with n 
vertices, the algorithm reports the correct answer in 0(n log n) time in the real- RAM model of computa- 
tion [18j . It constructs offsets with increasing radii in three stages; the intermediate shapes arising during 
the computation are in general more difficult to offset than polygons, as they are bounded by straight-line 
segments and "indented" circular arcs (namely, the shape is locally on the concave side of the arcs). The 
main observation is that for certain classes of such shapes, these circular arcs can be ignored when computing 
the next offset (see Thcorcm[5]for the precise statement). This observation bounds the time required by each 
offset computation by 0{n\ogn), which is the key to the efficiency of the decision algorithm. Our proof is 
constructive, that is, if a solution exists it can be computed with the same running time. 

The computation of the exact decision procedure requires the handling of algebraic coordinates of consid- 
erably high degree. As an alternative, we give an approximation scheme that works exclusively with rational 
numbers. The scheme proceeds by replacing the offset disks by polygonal shapes of similar diameter, whose 
precision is determined by another parameter 5 < e. We prove a bound A that depends on e, the minimal 
e for which the answer to the decision problem is YES, such that the rational approximation returns the 
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exact result for all 5 < A. If the input shape is found to be deconstructible, this algorithm also outputs a 
solution. The computation of e up to any desired precision is still possible. We believe that our investigation 
of the relation between S and e is of independent relevance, mostly to the study of certified algorithms that 
approximate geometric objects with algebraic coordinates by means of rational arithmetic. 

The deconstruction problem leads to natural optimization questions: if Q and r are given, how to compute 
e, the minimal tolerance for which a solution exists? Similarly, if Q and s are given, what are the possible 
radii such that a solution exists? For the first question, we provide a certified and efficient solution based on 
binary search, using the rational approximation algorithm. For the second question, we prove that the set of 
possible radii forms an interval and propose an algorithm to compute it. We also provide a heuristic to find 
a reasonable radius r if both r and e are unknown. 

For a convex shape Q with n vertices, we reduce the running time for solving Problem 1 to the optimal 
0{n) (in the real- RAM model). Moreover, we describe a greedy algorithm within the same time complexity 
that returns a solution P* which minimizes, up to one extra vertex, the number of vertices among all solutions, 
if there are any. Our algorithm technically resembles an approach for the different problem of finding a vertex- 
minimal polygon in the annulus of two nested polygons [T]. We also remark that the r-offset of P* has a 
tangent-continuous boundary and therefore constitutes a special case of an arc-spline approximation of Q 
where all circular arcs have the same radius. 

Organization. We describe an exact decision algorithm for the deconstruction problem (solving Prob- 
lem 1 above) in Section[21 In Section|3]we describe a rational-approximation algorithm for the deconstruction 
problem. Both algorithms output a solution in case the input is deconstructible (solving Problem 2). Sec- 
tion 2] discusses the optimization problems. For convex input. Section [5] exposes a specialized deconstruction 
algorithm and the computation of an almost vertex-minimal solution. We conclude in Section [5] by pointing 
out open problems. 



2 The Decision Algorithm 

For a set X C denote its boundary by dX and its complement by := M.'^ \ X. A polygonal region 
or polygonal shape X C is a set whose boundary consists of finitely many line segments with disjoint 
interiors. The endpoints of these straight-line segments are the vertices of the polygonal region. We assume 
henceforth that the input shapes that we deal with are bounded (but not necessarily connected). Although 
the techniques seem to go through also for unbounded shapes, this assumption simplifies the exposition and 
is sufficient for the real-life applications we have in mind. For two sets X and we denote their Minkowski 
sum hy X (BY := {x + y \ X & X , y £ Y}. With d{-, •) the Euclidean distance function, and any c e ]R^,r G R, 
we write Dr{c) := {p € | d{c,p) < r} for the (closed) r-disk around c, and Dr := Dr{0) for the r-disk 
centered at the origin. The r-offset of a set X, offset(A, r), is the Minkowski sum X Dr- 

For p S and X a closed set, we write d{p,X) := mm{d{p,x) \ x €E X}. The (symmetric) Hausdorff 
distance of two closed point sets X and Y is H{X, Y) :— niax{max{(i(x, Y) \ x G X}, max{d{y, X) \ y g Y}}. 
We say that X is e-close to Y (and Y to X) if H{X,Y) < e, which can also be expressed alternatively: 

Proposition 1. For X, Y closed, X is e-close to Y if and only ifYC offsct(A, e) and X C offsct(Y, e). 

Decision algorithm. We fix r > 0, £ > 0, and a polygonal region Q, and consider the following question: 
Is there a polygonal region P such that Q and the r-offset of P have Hausdorff-distance at most e? First of 
all, we can assume that r > e; otherwise, we can choose P :— Q, because oSset{Q,r) and Q have Hausdorff- 
distance at most e. We define another operation, r -inset (a.k.a. "erosion"), which is computationally similar 
to an offset: 

Definition 2. For r > 0, and X C M^, the r-inset of X is the set inset(A, ?■) := offset(A'-^, r)*^ = 
{xeR^l Drix) C X} . 

We are now ready to present the decision algorithm: 
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Algorithm 1 Decide(Q, r, e) 

(1) Q,^oSsct{Q,e) 

(2) n ^ inset((3e,r) 

(3) Q' ^ offset (n, r + e) 

(4) if Q C Q' then return YES else return NO 



We next prove that Decide (Algorithm [T]) correctly decides whether Q is e-close to some r-offset of a 
polygonal region. A first observation is that for any polygonal region P, offset(P, r) C if and only if P C n. 
This is an immediate consequence of the definition of the inset operation. This shows that for any offset (P, r) 
that is e-closc to Q, P must be inside 11. Moreover, it shows that any choice of P C n already satisfies one 
of Proposition [l]s inclusions. It is only left to check whether Q C offset(offset(P, r), e) — offset(P, r + e). We 
summarize: 

Proposition 3. Q is e-close to offset(P, r) if and only i/ P C 11 and Q C offset(P, r + e). 

To prove correctness of Decide, we have to show that Q C offset(n, r + e) already implies that there also 
exists a polygonal region P C n with Q C offset(P, r + e). The main difficulty in proving this is that 11 is not 
polygonal in general; we have to study its shape closer to prove that we can approximate it by a polygonal 
region, maintaining the property that the offset remains e-close to Q. 

The shape of offsets and insets. For a polygonal region Q, it is not hard to figure out the shape of 
Qe = oSset{Q,e): It is a closed set bounded by straight-line segments and by circular arcs, belonging to a 
circle of radius e. It is important to remark that all circular arcs are bulges: 

Definition 4. Let A C &e a closed set with some circular arc 7 on its boundary. Then, 7 is called a 
dent with respect to X , if each line segment connecting two distinct points on 7 is not fully contained in A. 
Otherwise, the arc is called a bulge. 

We call A a bulged (resp. an indented^ region with radius r, if dX consists of finitely many straight-line 
segments and bulges ( resp. dents ) that are all of radius r, interlinked at the vertices of the region. 

Note that a bulged region (left) is not necessarily convex. The r-offset of a polygonal 
region P is a bulged region with radius r. The heart of this section is Theorem[5]showing 
that the same also holds if P is an indented region (right) with radius smaller than r: 

Theorem 5. Let P be an indented region with radius ri, and let r2 > ri. Then, there 
is a polygonal region P^ C P such that offset(P, r2) — offset(PL, r2). In particular, 
offset(P, is a bulged region with radius r2. 




Figure 2.1: (a) The (extended) linear cap split by the polyline £, (b-c) the two cases in the proof of Theorem^ 
Proof. After possibly splitting circular arcs into at most four parts, we can assume that 

each circular arc spans at most a quarter of the circle. For such a circular arc, we define its endpoints by xi 
and X2, and denote the linear cap of the circular arc as the (closed) indented region enclosed by the circular 
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arc, and the two lines tangent to the circle through xi and X2 (the shaded area in Figure [?.lap . The extended 
linear cap is the (polygonal) region spanned by the two tangents just mentioned, and the two corresponding 
normals at xi and X2- Clearly, the normals meet in the center of the circle that defines the arc. 

We iteratively replace a indented arc of an indented region P' with radius ri (initially set to P) by a 
polyline £ ending in the endpoints of the circular arc, such that £ does neither leave P' nor the linear cap of the 
circular arc, and such that other boundary parts of P' are not intersected. This yields another indented region 
P" with radius ri , where one indented arc is replaced by a polyline, as depicted in Figure I2.1al Iterating 
this construction, starting with P, until all indented arcs are replaced, we obtain a polygonal region P^. 

We show that in each iteration, the r2-offsets of P' and P" are the same. For that we consider any point 
x' E P'\ P" , in the region that is cut off by P", and consider y ^ x' + v' for an arbitrary v' E 13^2 ■ We show 
that in all cases, y can also be written hy y — x" + w", with x" E P", and v" E D^.^. 

Since the circular arc spans at most a quarter of the circle, it is easily seen that D^^ {xi) U ZJ^i {X2) covers 
the whole extended linear cap. Therefore, for any y that lies within the extended linear cap, selecting x" = Xi 
or x" = a;2, we get y = x" + v" with v" E Dr^- 

We distinguish two other cases: for y that lies outside of the extended linear cap v' — x'y crosses either 
£ or the circular arc. In the former case, we can simply pick the crossing point as x" , and set v" E Dr^ 
accordingly (Figure I2.1bp . In the latter case, let us denote the crossing point as x* (Figure 12. Ic^ . We 
consider the set of points that is closer to x* than to xi and X2- Clearly, that region is bounded by the two 
corresponding bisectors, which meet in the center of the circle that defines the circular arc and is therefore 
completely contained within the extended linear cap. It follows that y is closer to one of the endpoints of 
the arc, say xi, than to x* . Selecting x" = xi we ensure y is closer to x" than to x' , which proves that 
y — x" + v" with some v" E Z?,^ in this case as well. □ 

The proof of Theorem [5] implies that offset (P,r2) for such a region P is completely determined by the 
offset of its linear segments, and the offset of the endpoints of circular arcs: the interior of the indented 
circular arcs can be ignored. 

Corollary 6. Algorithm]^ ('DEClDEj returns YES if and only if there exists a polygonal region P such that 
offset(P, r) is e-close to Q. 

Proof. Qe is a bulged region with radius e. Therefore, is an indented region with the same radius. 
Since r > e, Theorem [5] implies that oSset{Q^,r) is a bulged region with radius r, and so, offset(Q^, r)'" — 
inset((5e,r) = 11 is an indented region with the same radius. Using r + e > r and applying Theorem [S] once 
more, there exists a polygonal region P C n such that offset(n, r + e)= offset(P, r + e). It follows that, if the 
algorithm returns YES, there is indeed a polygonal region P whose r-offset is e-close to Q. If the algorithm 
returns NO, it is clear that no such polygonal region can exist. □ 

Theorem 7. Let P be an indented region with radius ri having n vertices, and assume r2 > ri. Then, 
offset(P, has 0{n) vertices and it can he computed in 0{n\ogn) time. 

Proof. By Theorem [SJ it suffices to consider a polygonally bounded Pl instead of P. We use trapezoidal 
decomposition of P to construct such a Pl with only 0{n) vertices. The Voronoi diagram of Pl's vertices 
and (open) edges can be computed in 0(n log n) time and has size 0{n) |22| . From it, the offset polygon 
with the same asymptotic complexity can be obtained in linear time |13j . □ 

Corollary 8. Algorithmic ('DEClDEj decides e-closeness with O(nlogn) operations. 

Proof. Apply Theorem [7] in each step of Algorithm [TJ The fourth step runs in 0(n log n) time as well using 
a simple sweep-line algorithm. □ 

Note that Hi,, if constructed for 11 as in the proof of Theorem [5] during step (3) of Algorithm [l] is a 
solution to the deconstruction problem if Decide returns YES. 
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3 Rational Approximation 



A direct realization of Algorithm [T] runs into difBculties since vertices of the resulting offsets are algebraic 
numbers whose degrees become high in cascaded offset computations. We next describe two approximation 
variants of Algorithm [Tl each producing a certified one-sided decision by approximating all disks in the 
algorithm with fc-gons. In order to make guaranteed statements about the exact e-approximability by r- 
offsets, we have to approximate the disks by a "working precision" S which is even smaller than e. Recall 
that Dr is the disk of radius r centered at the origin. For a,b E M., a < b define Da,b to be a polygon with 
rational vertices whose boundary lies in the annulus Df, \ Da- In the approximation algorithms, every disk is 
replaced with such a polygon lying inside a 5-width annulus. 

Interior approximation. In the first part of our algorithm, we ensure that the final approximation of 
Q' (see line (3) of Algorithm [T]), called Q' , will be a subset of the exact Q' . We achieve this by approximating 
Dg by Dg^s s when an offset is computed; and by approximating Dg by Dg g^^g when an inset is computed; 
see Algorithmic] 



Algorithm 2 ApproxDecideInterior {Q,r,e,S) 

(1) Q,^Q(SD^e with D, ^ D.^s^e 

(2) n<-((3£ (S Drj with Dr ^ Dr,r+S 

(3) Q' ^ n© Dr+e with Dr+e ^ Dr+e-S,r+e 

(4) if Q C Q', return YES, 
otherwise, return UNDECIDED 



Lemma 9. // ApproxDecideInterior {Q,r,e,S) returns YES, then Decide(Q, r, e) returns YES as well, 
which means that there exists a polygonal region P such that offset(P, r) is e-close to Q. In particular, P := 11 
is a solution to the deconstruction problem. 

Proof. Compare the execution of Algorithm [2] with the corresponding call of its exact version. Algorithm [TJ 
It is straight-forward to check that for any 6, Qe C Qg, 11 C 11, and Q' C Q' . The last inclusion shows that 
if Q C Q^, also QCQ'. □ 

Definition 10. For fixed Q and r, define e :— inf{e | Decide((5, r, e) returns YES}. 

Note that e E [0,^], and that Decide((3, r, e) returns YES for every e > e and returns NO for every 
e < e. We do not have a way to compute e exactly. However, we show next that ApproxDecideInterior 
(Q,r,s,d) returns YES for every e > e for S small enough, and that the required precision S is proportional 
to the distance of e to i. 

Theorem 11. Let e > i, and S < Then, ApproxDecideInterior {Q,r,e,S) returns YES. 

Proof. Let Eq be such that e < Sq < Sq + 26 < e. Let Qea, n and Q' denote the intermediate results of 
Decide {Q,r,eo) and let Q^, U, Q' denote the intermediate results of ApproxDecideInterior {Q,r,e,6). 
By the choice of Eq, YES is returned, and thus Q C Q' . The theorem follows from Q' C Q', which we prove 
in three substeps: 

(1) oSset{Q,„,6) C Ol: Indeed, 

oSset{Qeo , S) = offset((5, eo + S)C offset(Q, e - 5) c Q (BDI ^ QI- 

(2) n C 11: Starting with (1), we obtain 

offset(g,„,J) C Ql 
inset(offset((3eo,J),r + (5) C n . 
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We use the general fact inset(ofFset(yl, r), r) ^ A to obtain: 

inset(offset((3e(,, (5), r + (5) = inset(inset(offset((5eo i ^ inset((3eo, r) =11. 

(3) Q' C Qf; Using (2), we have that 

offset(n, r + e-5)=U(B Dr+eS C fi ® IvTe = Q' ■ 
Note that r + e — (5>r + eo, and therefore, offset(n, r + e — 5) D offset(n, r + eo) = Q' ■ D 

Exterior approximation. In Algorithm |3l we ensure that Q' becomes a superset of the exact Q' by 
appropriately choosing approximate disks. Specifically, we approximate Dg by Dg.s+s when an offset is 
computed, and Ds by Ds-s,s when an inset is computed. Not surprisingly, we get a certified answer in the 
other direction, and a certified answer is guaranteed when S is sufficiently small. The proofs of the following 
two statements are similar to Lemma [9] and Theorem [TT] and thus omitted. 

Algorithm 3 ApproxDecideExterior {Q,r,e,S) 
(I) Ql^Q(3lD, with D, ^ D,^,+s 

{2)Il^iQs QDr) with Dr ^ Dr~S,r 

(3) Q' ^ ft ® Dr+e with Dr+e ^ Dr+s,r+e+5 

(4) if Q C Q', return UNDECIDED, 
otherwise, return NO 



Lemma 12. // ApproxDecideExterior {Q,r,e,5) returns NO, then Decide(Q, r, e) returns NO as well, 
which means that there exists no polygonal region P such that offset(P, r) is e-close to Q. 

Theorem 13. Let e < e and S < Then, ApproxDecideInterior {Q,r,e,S) returns NO. 

In combination with Theorem [TTJ it follows that the exact answer can always be found toi S < A :— 
by combining ApproxDecideExterior and ApproxDecideInterior. We display the complete rational 
approximation algorithn:0 for later reference: 



Algorithm 4 ApproxDecide {Q,r,e,S) 

(1) if ApproxDecideInterior {Q,r,e,S) = YES, return YES 

(2) if ApproxDecideExterior {Q,r,e,S) = NO, return NO 

(3) Otherwise, return UNDECIDED 



Complexity analysis. The main task is to bound the number of vertices of 
Da.b- We will create a Da^b with the additional property that all vertices lie on 
dDb- As depicted on the right, two such points on Dt are connected by a chord 
of the boundary circle that does not intersect Da if and only if the angle induced 
by the two points is at most ijj := 2arccos|, or equivalently, the length of the 
chord is less than 1\/W~~a} . Note that we need at least ^ points on dD\, for a 



valid Da.h 



as easily shown by L'Hopital's rule. 




Rational points on dD\, can be constructed for an arbitrary t G Q as Qf := 
(^T^Iti ^i^It) [3 . For some positive z e Z, we define Pi := Qi/^ for i = 0, . . . , z. 

Lemma 14. For every i = 1, . . . , z — 1, the chord Pi^iPi is longer than the chord PiPi^i. In particular, the 
length of each chord is bounded by the length of PqPi which is shorter than —. 



^As in the Decide case for e = the return value is NO, and for e = r it is YES with P = Q. 
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Proof. W.l.o.g., we assume 6 = 1 for the proof, since the chord length 
scales proportionally when scaling the circle by a factor of h. The point 
Qt = (i^fi-, can be constructed geometrically as the intersection 

point of dDb with the line through S = (—1,0) and slope t (see the figure 
on the right). In particular, the line SPi has slope ^; we let Ti denote the 
intersection point of that line with the line a; = 1. We observe that the 
segment TiXi^i has length |, and that STi < ST^+i for i = 0, . . . , z — 1. 

We are showing next that the chord Pi-iPi is longer than PiPi+i. For 
that, we consider the triangle SDi-iDi+i, and its bisector at S. This bisector intersects the line a; = 1 at some 
point B. By the Angle Bisector theorem, B divides the segment Ti-iT^+i proportionally to the corresponding 
ST' BT' 

triangle sides, that is, /ir,^^^ — -nrrr—^- Because the left-hand side is smaller than 1, it follows that BTi^i 
bli+i Dli+i 

is shorter than ST^+i. Therefore, B lies below T^, and therefore, the angle ai^i — ZTj_iSTi = ZPi-iSPi is 
larger than = /.TiSTi+i = ZPiSPi+i. But the chord lengths Pi-iPi and PiPi+i are defined by 2sin(Q;i_i) 
and 2sin(Q;i), respectively, which proves that the chord lengths are indeed decreasing. 

Finally, by Thales' theorem, the triangle SP^Pi has a right angle at Pi. Therefore, the longest chord 
PqPi is shorter than the segment Tf)Ti, which has length |. □ 

Note that all P^'s lie in the first quadrant of the plane and that Pq := (6, 0) and P^ :— (0, b). Therefore, 
we can subdivide the other three quarters of the circle symmetrically such that the length of each chord 
is bounded by ^, using 4z vertices altogether. To compute a valid -Da, 6, it suffices to choose z such that 



that is z > 



_b_ _ 

b—a ^ V b—a fc+a y b'^—c 

asymptotically optimal choice. We summarize the result 

Lemma 15. For a <b, a polygonal region Da.b cls above with 0{^J'^^) (rational) points can be computed 
using 0(y^ arithmetic operations. 

The Minkowski sum of an arbitrary polygonal region with n vertices and a convex polygonal region with k 
vertices has complexity 0{kn) and it can be computed in 0(nfc log^(nfc)) operations by a simple divide-and- 
conquer approach, using a sweep line algorithm in the conquer step jl4j . Using generalized Voronoi diagrams 
where the distance is based on the convex summand of the Minkowski sum operation |15j . we obtain an 
improved algorithm, which requires only 0{kn\og{kn)) operations. In combination with Lemma 1151 this 
leads to the following complexity bound for the two approximation algorithms. 

Theorem 16. Algorithm ApproxDecide requires 



b2 



b2- 



We choose zq 



b 

b-a 



indeed, since < a < 6, we have that 



62 



As stated above, we need at least r2(w^^) points, so Zq is an 



arithmetic operations with rational numbers. 



We remark that the 0{nlogn) bound for Decide refers to operations with real numbers instead. 

We have implemented the algorithms ApproxDecideInterior and ApproxDecideExterior using 
exact rational arithmetic using the CgaiH packages for polygons [9], Minkowski sums [21] and Boolean set 
operations • We demonstrate the execution of our software on two examples in Figures 13.11 and 13.21 



4 Searching e and r 

So far, we have assumed that both r and e are given as input parameters, and we posed the question of 
deconstructing a polygon with respect to these parameters. We now investigate three variants where r 
and/or e are unknown. Specifically, we ask, for some input polygon Q: 

■^The Computational Geometry Algorithms, |www ■ cgal ■ org| 



8 



l-^^^= M (b) £ = 5 ■ <5 = ^ ■ (c) £ = I ■ r, 5 = I ■ e {d)e = ^-r,5=^ 



Figure 3.1: Dependency of the algorithm outcome on e and S: The input polygon (wheel) appears in bold 
line. It is colored according to its approximability with the given parameters: green for YES, red for NO 
and yellow for UNDECIDED. The inset polygon TP and its approximate (r + e)-ofFset Q' are drawn in green 



and cyan respectively. Their outer- approximation counterparts are drawn in red and magenta. Figures (a) 
and 1(b) I demonstrate how when e is tightened from ^ • r to ^ • r, with the same r and S, the decision result 



changes from YES to NO. The green polygon inside the input polygon in (a) is a possible r-offset solution. The 
magnification in |(b)| highlights the area of the input polygon that does not fit inside the outer ^-approximation 



(in magenta) of maximal possible (r -|- e)-offset. Figures (c) and |(d)| show how when S is decreased from j ■ e 



to ■ e, for the same r and e, the decision result changes from UNDECIDED to NO, namely in the latter case 
the algorithm is able to produce a certified negative answer. 




Figure 3.2: A map of Kazakhstan, represented as a polygon Q (in bold blue) with 1881 vertices, is ap- 
proximable for £ = 1 ' ^^"^ ~ I ' ^- ^ solution polygon P (in green) has 335 vertices. Offset(P,r) 
(shown as lightly-shaded gray r-strip around P) is inside the e-offset of the input Q by construction. The 6- 
approximation of the e-offset of Offset(P,r) (as computed in line (3) of ApproxDecideInterior {Q, r, e, S)) 
is drawn in cyan and has 261 vertices. Since the cyan polygon contains Q, the Offset(P,r) and Q have Haus- 
dorff distance of at most e, that is, Q is approximable and P is a solution. Approximability computation 
took 3.868 seconds in this case on a 3GHz Intel Dual Core processor. The magnification on the left highlights 
some cavities in the input polygon that have no effect on the Hausdorff distance within this tolerance e. The 
magnification on the right demonstrates a sharp end that would prevent Q's approximability with a tighter 

£. 
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1. Given r, what is e, the infimum of ah e- values such that the deconstruction problem has a sohition 
(compare Definition [TU)) ? 

2. Given e, what is the set of radii for which the deconstruction problem has a solution? 

3. Given neither r nor e, how to choose them in a "reasonable" way to obtain a solution? 

Whereas the first two questions are formally posed, the third one is of a rather heuristic nature. In all three 
cases, we also ask for computing some polygonal shape P that approximates the solution of the deconstruction 
problem for the given set of parameters. 

We discuss the posed questions in the remainder of this section. Our main tool will be the decision 
algorithm for fixed r and e as described earlier. Because we aim for a practical algorithm, we formulate our 
approach using the rational approximation algorithm from Section |3l We have implemented the proposed 
algorithms; the example at the end of this section has been produced with our implementation. 

Searching for e. If we use the exact decision procedure Decide, it is straight-forward to approximate e 
to arbitrary precision A employing binary search: Start with the interval [0,r] and choose e as the midpoint 
of the interval. If Decide((3, r, e) returns YES, recurse on the left subinterval, otherwise, on the right one. 
Obviously, the interval width is halved in every step, so 0(log(^)) steps are necessary. Let ea denote the 
£ approximation, s.t. — e < A. We demonstrate next that we can achieve the same approximation 
and produce with it a solution to the deconstruction problem using the rational approximation version 
ApproxDecide. 

Let |/| denote the width of / henceforth on and consider the pseudocode given in Algorithm^ It computes 
an interval / of width at most A that contains e. 



Algorithm 5 ApproxSearchEps(Q, r. A) 

(2) v^fhile |/| > A do 

(3) eno -i— left cndpoint of /, Sycs -s— right endpoint of / 

(4) e,„id^ ^^^4^' '^^ ^ 

(5) res ApproxDecide((5, r, Emid, 5) 

(6) if res — YES then / •(— [£„□, Emid] 

(7) otherwise, if res — NO, then I ^ [£mid,eycs] 

(8) otherwise, (res = UNDECIDED), then / ^ [emid — ^-j-,£mid + 

(9) end vifhile 

(10) return / 



We prove the invariant that e & I after each iteration of the while-loop, implying correctness of the whole 
algorithm. Trivially, e £ [0,r], and the invariant is obviously maintained if ApproxDecide((5, r, Emid) ^) 
returns YES or NO. For the case of UNDECIDED, recah that ApproxDecide is a combination of the two 
one-sided approximation algorithms ApproxDecideExterior and ApproxDecideInterior, and both 
returned UNDECIDED. Theorem [Til and Theorem [13] imply therefore that 

^ | £mid - £ I 

8 - ' 2 

It follows that £ g [£mid — ^ , Emid + ^] which provcs that the invariant is maintained also in this case. 

We next compute a solution P for the deconstruction problem for Q, r and £a- Recall that if Approx- 
Decide returns YES, the algorithm computes a solution for the deconstruction problem as a by-product. Let 
/ ApproxSearchEps(Q, A) be the approximation interval for e and let £ycs denote the right endpoint 
of /, that is £ < £yos. We call ApproxDecide((5, r, £yes, ^-^)- If the result is YES, then £a = £yos and the 
polygon computed by ApproxDecide is a solution. Otherwise let us choose £a = £yos + y and produce a 
solution by calling ApproxDecide(Q, r, £a, -f-). Since the result of ApproxDecide((5, r, £yos, ^) was UN- 
DECIDED we conclude from Theorem 1111 that e > Eyes — that is £a = £yes -|- ^ is indeed A-approximation 
of £. The cah to ApproxDecide((3, r, £a, -f-) is bound to yield YES because if it returned UNDECIDED, we 
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would have that e > ea ^ 2-j = eyes + -j, a contradiction to e < £yos- So, the polygon computed in this call 
is a solution. 

An overall complexity analysis of approximating i (and computing a solution) is relatively straightforward: 
/ is obviously halved in every iteration, so it takes 0(log(-^)) iterations to approximate e. Every iteration is 
bounded by the complexity given in Theorem 1161 We omit further details of the proof: 

Theorem 17. Approximating e to a precision A > requires 

vis r 6 

0{n—\ — ■ \og(n—\ — )) 
A V A A V A 

arithmetic operations with rational numbers. 

Searching vaUd radii. We assume now that Q and e are given, and discuss the question of what is the 
set R of radii such that the deconstruction problem has a solution. A priori, it is not clear what is the shape 
of R, but we will prove that it is an interval of the form [0, r*]. Having established this, we can apply another 
variant of binary search to approximate the extremal value r* . 

In order to prove that R is an interval, we prove first that the deconstruction problem can always be 
solved for Q, r and i. In other words, we can replace the infinium in Definition [10] by a minimum. The proof 
relies on two properties of infinite intersections of insets and offsets that we show first. 

Lemma 18. Let (yli)^^^ be a sequence of closed sets in M^. Then 

oo oo 

inset(P| Ai, r) = P| inset (A^, r). 

1=0 i=0 

Proof. The fact follows readily from the definition of insets: If a e inset(p|^Q Ai,r), then Dr{a) is contained 
in Pli^o particular, it is contained in Ai for every i which proves one inclusion. The other direction is 

similar. □ 



Lemma 19. Let (Ai)igN be a sequence of closed sets in with Aq ^ Ai ^ . . .. Moreover, let (Ai)igN be a 
monotonously decreasing sequence of real numbers that converges to A G M. Then 

oo oo 

offset(P| A,,A) = Pi offset (^i. A,). 

i=0 i=0 

Proof. The "C" inclusion is straight forward, so we concentrate on the "3" part. Fix some b ^ f] offset(^i, Ai). 
For every i G N, there exists some G Ai such that [b—ai) G D\^. Now, the sequence (6— ai)igN is a bounded 
sequence in M? (bounded by -Daq) ^J^d therefore has a convergent sub-sequence by the well-known Bolzano- 
Weierstrass Theorem. Let r denote the limit point of this subsequence. In particular, the corresponding 
subsequence of (ai)igN converges to a := 6 — r. We show that a G a-^^ '"' ^ which suffices to prove 

the claim. 

Assume that a ^ H ^i- Then, there is some no such that a ^ An^. Since A„„ is closed, c?(a, An^) =: e > 0, 
where d is the Euclidean distance function. Moreover, because each An with n > ng is included in An„, 
d{a, An) > e. Because {ai)i^ti converges to a, we can find some N > uq such that d{a,aN) < e. However, 
Cat G Am, so 

d{a, Am) < d{a, om) < e = d{a, An^) < d{a, Am), 
which is a contradiction. The fact that r G D\ follows by a similar argument. □ 

Theorem 20. For arbitrary Q and r, and e as from Definitional^ there exists a solution to the deconstruction 
problem. 

Proof. Because of Corollary |6l we need to prove that 

Q Q Q'i '■— offset(inset(offset(Q, e), r), r -I- e). 
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Let (£i)igN be a monotone decreasing sequence of real numbers that converges to e. Because > £ for each 
i, Decide return YES for each Ci, which is equivalent to 

Q Q Q'i ■= offset(inset(offset(Q, e^), r), r + Si). 

It is therefore sufRcient to prove that 

oo 

1=0 

For that, we apply Lemma fT9l on the constant sequence ((5)igN and on (ei)^^^ to obtain 

oo 

offset(Q,£) = P) offset((5, £i). 

Applying Lemma [TH] yields 

CxD 

inset (offset (Q, e), r) — (inset(offset((3, e^), r), 

1=0 

and applying Lemma [TOl for the sequences (inset(offset((3, e^), and (r + £i)ieN yields 



offset(inset(offset(Q, e), r),r + e) = \ \ offset(inset(offset((3, e^), r), r + Si) . 




□ 

Since r is no longer fixed, we now consider e as a function from to M'*' depending on r; abusing 
notation, we will let e by itself denote this function from now on. 

Theorem 21. e is a monotone increasing function. Moreover, e is Lipschitz-continuous with Lipschitz factor 
1. 

Proof. We prove monotonicity first: Let a < b denote two radii and Ea = s{a), £b = We will show that 
Decide (a,eb) returns YES, which proves that Ea < £&■ 

We compare the intermediate results of Decide (a, e&), denoted by Si, to those of Decide (6, e;,), denoted 
by B,: 

5*1 = offset((5, £f,) Bi — offset(Q, £(,) = 5*1 

5*2 = inset(S'i, a) i?2 = inset(i?i, h) = inset (inset (5*1, a),b — a) = inset(S'2, b — a) 

5*3 = offset(S'2, a + Eb) — offset(i32, & + = offset (offset (i?2, & — a), a + Eb) 

Since 

ofFset(i32, b — a) = offset (inset (^2, b ~ a) ,b ~ a) Q S2 

it follows that B^ C 53. Because Decide (6,eb) returns YES by definition, it holds that Q C ,63, therefore 
Q C S'3 and Decide {a,Eb) also returns YES. 

For Lipschitz continuity, let a < 6 be such that 5 — a < (5, and again Ea = e(a), £{, — £{b). We show that 
£b — £a < <5. There exists a polygonal region P that is a solution to the deconstruction problem for Q, a and 
Ea- In other words, 

i/(offset(P,a),Q) < Ea- 

Because of the general fact 

H{A, B)<E^ iJ(offset(A, 6),B) <E + 6, 

and since b < a + S wc have that 

F(ofrset(P, 6), Q) < H{oSset{P, a + 6), Q) = (offset (offset (P, a), S),Q)<Ea + 5. 
Therefore, P is a solution for the deconstruction problem for Q, b and Ea + S, so Eb < Ea + S. □ 
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It follows from the monotonicity and Theorem [50] that R is an interval which has as its left endpoint. 
Thus, computing R reduces to finding the maximal r* > such that e(r*) = e. 

Using the exact decision procedure, we can perform a binary search similar to that for approximating 
e: First, we compute an interval [0,r] containing r* . Since Q is finite, we can take r to be the radius of 
the smallest enclosing circle of Q plus e. Then, we start the iterative process, deciding on the left or right 
subinterval depending on the result of Decide for Q, e and the midpoint of the interval. 

What if we are using ApproxDecide instead? Unlike Algorithm^ we can no longer guarantee that every 
execution of the approximation algorithm halves the search interval, because a return value UNDECIDED does 
not bound the distance of the current radius r to the critical value r*. Instead, we propose the following 
scheme: For an interval / with midpoint r, ApproxDecide is called with some 5, initially set to |. If it 
returns UNDECIDED, S is divided by 2 and ApproxDecide is recalled. Eventually, the algorithm returns 
YES or NO, and the interval / can be halved. 

Let R' be the preimage of e under e. Note that R' is an interval (which may consist of only one point). 
The algorithm from above is guaranteed to converge to some r € R'. However, if R' contains more than 
one point, it is not guaranteed to converge to the maximal one (because it gets stuck in an infinite loop as 
soon as the query value r lies in R'). One way of avoiding this infinite loops is to decrease S only to some 
threshold and choosing another query value r from the interval if no decision was made. Nevertheless, we 
have not found an algorithm with the formal guarantee of converging to the largest value in R' eventually. 

Searching for both r and e. We finally consider the question of how we can find a reasonable 
choice of r, e and a polygonal region P, such that P is a solution for the deconstruction problem for Q, r, 
and e. The meaning of "reasonable" depends on the application context, and possible prior information (for 
instance, a range of possible offset radii). We offer a basic generic approach and justify our choice with an 
example. 

Generally, we expect from a reasonable pair (r, e) that e is small. So, in order to judge whether a good 
solution exists for radius r, we consider e(r). However, e(r) being small (or equivalently, being large) 
is not a good criterion, because i is monotone increasing according to Theorem I2H so r = would always 
be the best solutionis In order to remove the bias towards small radii, we scale the objective function and 
consider 

J : M+ ^ M+,r 

e{r) 

Note that J is well-defined on the positive axis and continuous. Moreover, we can approximate the graph of 
J in any finite interval of r- values by choosing a sample of the interval and approximating i at each sample 
value using Algorithm [5] 

We demonstrate by an example that the local maxima of J yield radii that lead to good deconstruction 
results. Consider the polygonal region defined in Figure B.lai and its (approximated) J-graph in Figure l4.1bl 
We can identify two local maxima ri and ; we have plotted the corresponding solutions in Figure 14.21 
Indeed, we see that for the large radius ra, we obtain a relatively simple solution whose offset blurs away the 
spikes of Q. For the smaller local maximum at ri, we obtain a solution with more details such that the spikes 
can be approximated almost perfectly with the given radius. In contrast, the shape at the local minimum r2 
combines the disadvantages of the two discussed cases: the solution is similarly complicated as the ri-solution 
(it contains flattened versions of all the spikes), but its approximation quality is not significantly better than 
for the rs solution which achieves the same with a much larger radius. 

5 Deconstructing Convex Polygons 

Assume that the input Q to Algorithm [T] is a convex polygon. We first improve the decision algorithm such 
that it runs in linear time (Algorithm IB]). Then we look for a polygon P with a minimal number of vertices 
(OPT) such that Q is e-close to offset(P, r). We give a simple linear-time algorithm that produces a polygon 
with at most OPT + 1 vertices. 

Lemma 22. If Q is a convex polygonal region, then II, as computed by Decide (Algorithm]^, is also a 
convex polygon, and it can be computed in 0{n) time. 

•^Note that this is also formally correct, because P ;= Q is the perfect solution for r = 0. 
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(a) The polygon (b) The J-graph approximation with A = 

Figure 4.1: The Flower polygon example: 10 samples per radius unit for r < 5. 




(a) J maximum at ri = 0.6 (b) J minimum at r2 = 1.4 (c) J maximum at = AA 




(d) e^iri) Si 0.082 • n = 0.049 ... (e) s^(r2) » 0.202 ■ r2 = 0.283 ... (f) EACfs) ~ 0.099 ■ ra = 0.439 . . . 



Figure 4.2: Flower polygon approximations for the (r, £A(r)) values at J-graph extrenia. In the upper row, 
approximations of the solutions are shown in green, and their r-offsets (eA(r)-close to the input) in dark 
blue. In the lower row the input polygon is shown in blue. The SA-width cyan stripe around the r-offset 
demonstrates the quality of the approximation. 
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Proof. Q is the intersection of the half-planes bounded by lines that support the polygon edges. Observe 
that n can be directly constructed from Q by shifting each such line by r — e inside the polygon, which shows 
that n is convex. For the time complexity, we divide the shifted edges of Q into those bounding Q from 
above, and those bounding Q from below (we assume w.l.o.g. that no edge is vertical). Consider the former 
edges; the lines supporting those edges have slopes that are monotonously decreasing when traversing the 
edges from left to right. We have to compute their lower envelope; for that, we dualize by mapping y — mx + c 
to (m, — c), which preserves above/below relations, and compute the upper hull of the dualized points. Since 
we already know the order of the points in their x-coordinate, this can be done in linear time using Graham's 
scan [2 [To]. The same holds for the edges bounding Q from below, taking the upper envelope/lower hull. □ 

Decide first computes 11 and checks whether Q C offset(n,r + e). We replace the latter step for convex 
polygons: Let qi, . . . ,qn be the vertices of Q (in counterclockwise order) and define Ki ~ Dr+e{qi), namely 
the disk of radius r + e centered at qi. We check whether all these disks intersect 11: 



Algorithm 6 ConvexDecide(Q, r, e) 

(1) Qe ^ offset(g,e) 

(2) n ^ inset((5e,r) 

(3) if if , n n ^ for all i = 1, . . . , n, return YES, 
otherwise return NO 



Lemma 23. ConvexDecide agrees with Decide on convex input polygons Q and runs in 0{n) time. 

Proof. For correctness, it suffices to prove that offset(n, r) is e-close to Q if and only if each Ki intersects 11: 
Indeed, if any Ki does not intersect 11, then has distance more than r + e to 11, so Q is not e-close to the 
offset. Otherwise, if each disk Ki intersects 11, offset(n, r -\- e) contains each vertex of Q. Since it is a convex 
set (as the Minkowski sum of two convex sets), it also covers each edge of Q. Thus, Q C ofFset(n, r + e), 
which ensures that Q is e-close to the offset by Proposition [3] 

For the complexity. Lemma [22| shows that the computation of H runs in linear time. We still have to 
demonstrate that the last step of the algorithm (checking for non-empty intersections) also takes a linear 
time. Let ei, . . . , Cm be the edges of 11 (with m < n). To check for an intersection of Ki with H, we traverse 
the edges and check for an intersection, returning NO if no such edge is found. However, if such an edge, say 
Cj was found, we start the search for an intersection of the next disk if^+i at e^, again traversing the edges 
in counterclockwise order. Using this strategy, and noting that Ki, . . . , Kn are arranged in counterclockwise 
order around H, it can be easily seen that we iterate at most twice through the edges of 11. □ 

Reducing the number of vertices. We assume that offset(n,r) is e-close to Q. We prefer a simple- 
looking approximation of Q, thus wc seek a polygon P C H whose offset is e-close to Q, but with fewer 
vertices than H. Any such P intersects each of the bulged regions of radius r + e: := A'^ n H, i = 1, . . . , n. 
We call these bulged regions H's eyelets. The converse is also true: Any convex polygon P C H that intersects 
all eyelets ki, . . . , k„ has an r-offset that is e-close to Q. 

The following observation is a simple consequence of Proposition [3l 

Proposition 24. // offset(P, r) is e-close to Q, and P C P' C then offset(P',r) is e-close to Q. 

We call a polygonal region P (vertex-)minimal, if its r-offset is e-close to Q, and there exists no other 
such region with fewer vertices. Necessarily, a minimal P must be convex - otherwise, its convex hull CH(P) 
has fewer vertices and it can be seen by Proposition [24l that ofi'set(CH(P), r) is also e-close to Q. By the 
next lemma, we can restrict our search to polygons with vertices on OH. 

Lemma 25. There exists a minimal polygonal region P C H the vertices of which are all on dH. 

Proof. We pull each vertex pi ^ dH in the direction of the ray emanating from 
Pi-i towards pi until it intersects dll in the point p'i (dragging pi's incident edges 
along with it); see the enclosed illustration. For P' — {pi, . . . ,pi-i,p'i,pi+i, . . . ,Pm): 
P C P' C n, offset(P', r) is e-close to Q by Proposition [Ml □ 

We call a polygonal region P good, if P C n, all vertices of P lie on dH, and P 
intersects each eyelet ki, . . . , k„. Note that any good P is convex. 
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Definition 26. For two points u, u' G 911, we denote by [u, u'] C 911 all points that are met when travelling 
along dll from u to u' in counterclockwise order. Likewise, we define half-open and open intervals [u,u'), 
{u, u'], {u, u'). 



Let Ki = Ki {^H he qi's eyelet as before. Consider Ki D dH. The portion of tlrat 
intersection set that is visible from qi (considering IT as an obstacle) defines a (ccw- 
oriented) interval [vi,Wi] C 911. We call Vi the spot of the eyelet ki. Finally, for u,u' G 
9n, we say that the segment uu' is good, if for all spots Vi € {u,u'), uu' intersects the 
corresponding eyelet k^. 




The figure above illustrates these definitions: The segment pp' is good, whereas pp" is 
not good, because V2 G {p,p"), but the segment does not intersect k.2- 



Theorem 27. Let P be a convex polygonal region with all its vertices on 911. Then, P 
is good if and only if all its bounding edges are good. 



Proof. We first prove that if all the edges of P are good, then P is good. It suffices to " ~ - 

argue that it intersects all eyelets ki, . . . , k„. Let pi, . . . ,pk be the vertices of P in counterclockwise order. 
Any spot Vi of an eyelet Hi either corresponds to some vertex pi of P, or lies inside some interval {pi,pi+i). 
Since pipi+i is good, it intersects k^. 

For the converse, assume that pipi+i is not good, which encloses with the interval {pi,pi+i) a polygonal 
region i? C n \ P. Hence, there is a spot Vi d R such that pipi+i does not intersect the eyelet Ki. It follows 
that the entire Ki is inside R (see the illustration above, considering pp" and K2). Thus, P H Ki = and so 
P cannot be good. □ 

For u e 9n, we define its horizon /i„ S 911 as the maximal point in counterclockwise direction such that 
that segment uhu is good. Consider again the figure above: The segment php is tangential to K2, so if going 
any further than hp on 911 from p, the segment would miss K2 and thus become non-good. 

Lemma 28. Let P be a good polygonal region, and u G 911. Then, P has a vertex p £ (u, hu]. 

Proof. Assume to the contrary that P has no such vertex, and let pi, . . . ,pi be its vertices on 911. Let pj be 
the vertex of P such that u € (pj,pj+i). Then, also hu € {pj,pj+i), because otherwise, Pj+i G (u, hu]- Since 
P is good, the segment pJpj+T is good, too. It is not hard to see that, consequently, both pju and wpj+i are 
good. However, the latter contradicts the maximality of the horizon hu- □ 

For an arbitrary initial vertex s G 9H, we finally specify a polygonal region P'* by iteratively defining its 
vertices. Set pi := s. For any j > 1, if the segment pjs, which would close P'^, is good, stop. Otherwise, set 
Pj+i := hp. . Informally, we always jump to the next horizon until we can reach s again without missing any 
of the eyelets. By construction, all segments of P* are good, so P* itself is good. The (almost-)optimality of 
this construction mainly follows from Lemma 1281 

Theorem 29. Let P be a minimal polygonal region for Q, having OPT vertices. Then, for any s G 9H, P* 
has at most OPT + 1 vertices. 

Proof. We first prove that P" has the minimal number of vertices among all good polygonal regions that 
have s as a vertex. Let s := pi, . . . ,pm be the vertices of P*. There are m — 1 segments of the form pihp^. 
By Lemma UHl any good polygonal region has a vertex inside each of the intervals {pi,hp^]. Together with 
the vertex at s, this yields at least m vertices, thus P'* is indeed minimal among these polygonal regions. 

Next, consider any minimal polygonal region P*. We can assume that all its vertices are on 9H by 
Lemma [25l If s is not a vertex of P* , we add it to the vertex set and obtain a polygonal region P' with at 
most OPT + 1 vertices that has s as a vertex. P" has at most as many vertices as P', so m < OPT + 1. □ 

As each visit of an eyelet requires constant time, the construction of a horizon is proportional to the 
number of visited eyelets, and there are only linearly many eyelets. Thus, we can state: 

Theorem 30. For an arbitrary initial vertex s, computing P* requires 0{n) time. 
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Proof. We prove that computing the horizon of a point u takes a number of operations proportional to the 
number of eyelets that are visited by the segment uhu- Let us consider an arbitrary u G 511. By rotating 
appropriately, we can assume, without loss of generality, that u lies on a vertical edge of 11 (or, if w is a 
vertex, that the next edge in counterclockwise; ordc;r is vertical), and that the edge is traversed top-down. 
The horizon is determined by the slope of the edge at u. Note that for each eyelet ki, . . . , there is an 
interval of slopes , • • • , In'^ such that the segment from u with slope A intersects k, if and only if A € . 
Furthermore, each single I^^^ can be computed with a constant number of arithmetic operations. Assuming 
that the next eyelet to be travelled from the current pi is Kj. we can iteratively compute the intersections 
Ij n /j+i n Ij+2, ■ ■ ■ until /j n . . . n Ij-^-k is empty. In this case, wc choose Ai := max(/j n . . . n Ij+k-i) as the 
slope for the next segment, which must be Pihp^ since it is good by construction, and any larger slope would 
produce a non-good segment. Based on this property, it is easy to show that computing needs a number 
of operations which is proportional to n, the number of eyelets. □ 

6 Open Problems 

We have shown how to decide whether a given arbitrary polygonal shape Q is composable as the Minkowski 
sum of another polygonal region and a disk of radius r, up to some tolerance e. Many related questions 
remain open, (i) Dcconstruction of Minkowski sums seems more difficult when both summands are more 
complicated than a disk; many practical scenarios may raise this general dcconstruction problem, (ii) It would 
be interesting to analyze the dcconstruction not only under the Hausdorff distance but for other similarity 
measures, such as the Frcchct or the symmetric distance, (iii) Can one remove the extra vertex when seeking 
an optimal (vertex minimal) polygonal summand P in the convex case, (iv) Finding an optimal or near- 
optimal polygonal summand in the non-convex case seems challenging, (v) As in polygonal simplification, 
we could also search for the polygonal region with a given number of vertices whose r-offset minimizes the 
(Hausdorff) distance to the given shape, (vi) The offset-deconstruction problem can be reformulated in higher 
dimensions. We consider especially the three-dimensional case to be of practical relevance. 
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